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We define a class of insulators with gapless surface states protected from localization due to 
the statistical properties of a disordered ensemble, namely due to the ensemble's invariance under 
a certain symmetry. We show that these insulators are topological, and are protected by a Z2 
invariant. Finally, we prove that every topological insulator gives rise to an infinite number of 
classes of statistical topological insulators in higher dimensions. Our conclusions are confirmed by 
numerical simulations. 
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Introduction. One common definition of a topological 
insulator is that it is a bulk insulator which has a gapless 
surface Hamiltonian that cannot be continuously trans- 
formed into a gapped one [TJ [2] ■ The surface states of a 
topological insulator are protected from Anderson local- 
ization and, since there is an anomaly associated with the 
bulk field theory of such a material [3H5] , they are also ro- 
bust against interactions as long as the latter respect the 
discrete symmetry of the system. Other possible descrip- 
tions of topological insulators arise naturally from classi- 
fications of topological terms in nonlinear sigma-models 
[BJ, from K-theory [7J, and even from string theory [5]. 

There are, however, several known examples of dis- 
ordered systems which have a surface Hamiltonian that 
can be continuously deformed into a gapped one, and 
yet protected against Anderson localization. Remark- 
ably, these systems are closely related to conventional 
topological insulators. One such example is a so-called 
weak topological insulator, a three-dimensional material 
made by stacking many layers of a two-dimensional topo- 
logical insulator. Its surface has two Dirac cones which 
can be coupled by a mass term, producing a gapped sys- 
tem. Despite that, Ringel et al. [5] have argued that 
since an odd number of weak topological insulator layers 
is conducting, its surface must always be metallic. This 
prediction was tested numerically [TO] and later explained 
by Fu and Kane [TT] in terms of Z 2 vortex fugacity of a 
corresponding field theory. Another example of a system 
avoiding localization is a topological insulator subject to 
a random magnetic field as long as the magnetic field is 
zero on average |12j . A random magnetic field opens a 
random sign gap in the surface dispersion, which drives 
the surface to a critical point of the Chalker-Coddington 
network model [13 . 

These two examples share one common trait: in or- 
der for the surface to avoid localization, the disordered 
ensemble has to be invariant under a certain symmetry: 
translation by a single unit cell in case of a weak topolog- 
ical insulator or time-reversal for a strong topological in- 
sulator with a random magnetic field. Another common 
feature is that the surface theories of both of these sys- 



tems are in the middle of a topological phase transition. 
We show that these two features define a broad class of 
systems, which we call statistical topological insulators 
(STI). An STI is an ensemble of disordered systems in- 
variant under a certain average symmetry, in which each 
system may have discrete symmetries. The elements of 
this ensemble have surfaces protected from localization 
due to the combined presence of the average symmetry 
and the discrete ones, if any. For example, in the case 
of a weak topological insulator the average symmetry is 
translation by one unit cell, while the discrete symmetry 
is time-reversal. 

Some STIs without disorder become topological crys- 
talline insulators, introduced by Liang Fu [THUS], or in 
other words, they have a gapless surface dispersion pro- 
tected by their crystalline symmetry. Nevertheless, we 
show that not all the topological crystalline insulators 
become STIs once disorder is added, and that the en- 
semble symmetry need not be crystalline, as in the case 
of a topological insulator in a random magnetic field. We 
also show that STIs are a true bulk phase: in order for the 
surface to become localized without breaking the symme- 
tries, the bulk must undergo a phase transition. This sig- 
nificantly extends the list of possible topological phases: 
since the bulk transition of an STI is a topological phase 
transition by itself, it is possible to construct a higher 
dimensional system with its surface pinned to the middle 
of an STI phase transition. Such a construction makes 
every single topological phase transition [6j [7J give rise to 
infinitely many higher-dimensional descendant topologi- 
cal phases. 

In the following, we start by defining the STI topologi- 
cal invariant when the average symmetry group that pro- 
tects it is Z2. Then we show how to build a tight-binding 
model for an STI in any dimension and symmetry class. 
Finally we illustrate the properties of STIs numerically 
in three 2D examples. 

Construction of an STI topological invariant. To de- 
termine the necessary conditions required to obtain a sta- 
tistical topological insulator, let us consider an ensemble 
of d-dimensional systems which are semi-infinite along 
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a) S = -1 b) Q = +1 c) AQ odd 




FIG. 1. (a), (b): Examples of possible patterns of surface 
gap closing as a function of the parameter A, interpolating 
between a system and its symmetry reverse, and the surface 
momenta k. Each gap closing contributes ±1 to the total 
change of the topological invariant. A nontrivial STI (a) has 
an odd number of gap closings at A = 0, while a trivial one (b) 
has this number even, such that the total change of topologi- 
cal invariant from X — 1 to A = — 1 may be 0. (c): the pattern 
of gap closings of an STI with two independent ensemble sym- 
metries. Since there is an odd number of gap closings that 
map on themselves under both symmetries, the change of a 
d — 1-dimensional STI invariant Q has to be odd. 



at least one spatial direction, such that they have one 
(d — l)-dimensional surface. We require that the Hamil- 
tonian Hi of each element of the ensemble is local and 
the correlation function of the matrix elements of this 
Hamiltonian is sufficiently short range. Additionally, we 
require that the bulk is insulating, while the surface may 
have surface states. We consider d > 2, so that both 
the surface and the bulk are self-averaging [16]. Let ev- 
ery ensemble element Hi be equally likely to appear as 
another one, zLUHiU^ 1 , obtained from Hi by applying 
a unitary or an anti- unitary Z 2 symmetry IA. Examples 
of such a symmetry are mirror, time-reversal, inversion, 
particle-hole, and chiral symmetries. 

Every element of the ensemble has the same discrete 
symmetries, allowing its surface to have a given topologi- 
cal invariant Qd—i- We consider elements of an ensemble 
of finite but large systems of size L closed with twisted 
periodic boundary conditions in all directions apart from 
the one perpendicular to the surface. In this way we may 
study d— 1-dimensional bands formed on the surface. In 
d > 2, the surface is macroscopic, so an insulating phase 
would have a definite value of Qd-i in the limit L — > 00. 

We may calculate the amount by which Qd-i changes 
under application of U. Since the d-dimensional bulk is 
gapped, this change equals the sum of the changes of 
Qd-i with each closing of the surface gap through a con- 
tinuous path which connects Hi and zLUHiU -1 . In order 
to make use of the average symmetry we choose this path 
H(X) such that H(l) = Hi and H(X) = ±UH{-\)U- 1 . 
The contribution of each gap closing to the topological 
invariant can be evaluated using the Green's function for- 
malism [18-20 , which allows to evaluate the topological 
invariant change Qd-i from the surface Green's function 
as long as the bulk stays gapped. Each gap closing con- 
tributes a change of ±1 to the topological invariant. 



All the gap closings belong to one of two groups: they 
either map on themselves under U, or on another gap 
closing, as shown in Fig. [lj,). The total change of the 
topological invariant can be expressed as a sum over these 
two types of gap closings, 

AQd-i = Qd-i(A = +1) - Q d -i(A = -1) 

= 2 J2 1n+ J2 (1) 

n:A„<0 n:A„=0 

where A n is the value of A where n-th gap closing is lo- 
cated and q n = ±1 is its contribution to AQ^-i. The 
parity of AQ d _i equals to the parity of the last term of 
Eq. ([I]). This parity should be the same for any element 
of the ensemble, since a path Hj — > Hi — > ±UHiU~ 1 — > 
zkUHjU^ 1 has the same number of gap closings at the 
symmetric point as the path Hi — > ^UHiU -1 . Further- 
more, this parity cannot be changed by perturbations 
that do not close the bulk gap, since the gap closings 
carry a finite topological invariant, and cannot be re- 
moved by local perturbations. A single gap closing can 
also not be shifted away from A = since this would 
break the symmetry. 

Let us consider the situation when AQ^-i is odd. This 
means that Hi and zLUHiU^ 1 may not have the same 
topological invariant. Since these are equally probable 
elements of the ensemble and Qd-i is a self-averaging 
quantity in insulating phases, this means that the surface 
may not be in an insulating phase. On the other hand, 
in the case when AQd-i is even, as shown in Fig. [IJd), 
the ensemble symmetry does not prevent the trivial in- 
sulating phase with Q = from appearing, and does not 
protect the surface states from localization. We thus de- 
fine Q = (— l) A< 2rf-i as the Z 2 topological invariant of an 
STI. The STI topological invariant Q is a bulk property, 
e.g. the parity of the mirror Chern number, the weak or 
strong topological insulator invariant. Nevertheless, an 
explicit evaluation of Q for large disorder strength is in 
general a hard problem, since it is necessary to find a 
Hamiltonian with exact symmetry U, which can be con- 
nected to an element of the ensemble without closing the 
bulk gap. 

The same construction works if we use an STI topo- 
logical invariant Q as the surface topological invariant 
Qd-i- We now consider a system with two indepen- 
dent ensemble Z 2 symmetries U\, lA^. Let Q be the par- 
ity of the changes of a topological invariant of a d — 2- 
dimensional edge of the surface under application of IA\ . 
We will now calculate the change of Q under applying a 
different symmetry Ui- Now, for each element of the en- 
semble Hi we need to consider a family of Hamiltonians 
#i(Ai,A 2 ), such that fli(Ai,A 2 ) = ±U 1 H{-X 1 , A 2 )Wf 1 
and Hi(\i,\ 2 ) — ±W 2 -ff(Ai, — A 2 )MT 1 . Using dimen- 
sional reduction arguments of Ref. [20], we obtain that 
the parity of AQ equals the parity of the total topologi- 
cal invariant carried by the gap closings inside the square 
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I Ai|, | Aa] < 1. Consequently, Q should change under ap- 
plying IA2 if and only if the number of gap closings that 
are mapped on themselves under application of U\ as well 
as U2 is odd, as shown in Fig. [TJ;). Repeating the same 
construction allows to construct STIs in dimension d us- 
ing an ensemble ZJ symmetry and a d — n topological 
invariant. 

STIs with mirror symmetry. Ensemble mirror sym- 
metry allows to construct STIs using any d-dimensional 
topological invariant Q. Consider an infinite set of d- 
dimensional layers stacked in the d + 1-st dimension and 
weakly coupled to each other. We alternate layers of A- 
and £?-type characterized by bulk topological invariants 
Qa and Qb respectively. We take Qb to be the inverse 
group element of Qa (i-e. Qb = ~Qa for Z invariants 
and Qb = Qa for Z2), so that coupling two such layers 
results in a topologically trivial system. Such a model has 
a reflection symmetry with respect to each layer. Once an 
onsite disorder with no long range correlations is added, 
this symmetry is broken, and becomes an ensemble sym- 
metry instead. 

For a sufficiently weak disorder we may choose the 
clean system as tt(\ — 0). Let us now evaluate the to- 
tal number of the fermi surfaces at the symmetric point. 
The general form of the Bloch wave function for this sys- 
tem is *(x) = exp(zfcr) (^ A , with ^ A and * B the 
wave function on A and B layers respectively, x the num- 
ber of the unit cell, and k Bloch momentum. At k = n 
the eigenstates of the Hamiltonian should also be eigen- 
states of mirror symmetry. These are exp(i7rx)(\I' J 4, 0) 
and exp(i7rx)(0, ^b), so they belong exclusively to ei- 
ther A- or £?-type layers. They also must carry a total 
topological invariant Qa and Qb at zero excitation en- 
ergy, since the layers are characterized by a non-zero bulk 
topological invariant. At a finite value of A the eigen- 
states of mirror symmetry couple, such that according 
to the dimensional reduction arguments of Ref. [20] the 
total topological invariant of the gap closings that map 
on themselves under the symmetry is Qa- Hence if each 
layer carries an odd number of topologically protected 
edge states, the layered system is a nontrivial STI and 
has Q = -1. 

Numerical Simulations. As an example of a two- 
dimensional STI obtained via a procedure similar to the 
one described above, we consider a stack of coupled Ki- 
taev chains oriented in the y-direction. The two dimen- 
sional lattice Hamiltonian of this system has the form 

H = (2t cos k y — V)cr z + Act y sin k y + aa x sin k x , (2) 

where o~i are Pauli matrices in Nambu space, and x, y are 
integer coordinates perpendicular and along the chain 
direction. The Hamiltonian parameters are chosen as 
follows: the normal hopping is t = 1, the onsite dis- 
order potential V is uniformly distributed in the inter- 
val [— e/2,e/2], the p-wave pairing strength within each 
chain is A = 1, and the interlayer coupling strength is 
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FIG. 2. Transmission through a stack of 49 coupled Kitaev 
chains, each chain 150 sites long with hard wall (green) and 
periodic (red) boundary conditions as a function of disorder 
strength e. Each point is obtained by averaging over a number 
of approximately 6 x 10 3 different disorder realizations, with 
error bars smaller than symbol sizes. The inset shows the 
sample geometry with Kitaev chains drawn as black lines, 
and with leads attached shown as gray rectangles. 



a = 0.45. This Hamiltonian belongs to symmetry class 
D with particle-hole symmetry H = —a x H*a x . It is also 
related to a Hamiltonian with an ensemble mirror sym- 
metry by a gauge transformation tp(x) — > (— 1) L^/ 2 J ip(x). 

The transmission in the x-direction of a stack of such 
Kitaev chains with periodic and hard wall boundary con- 
ditions in y-direction is shown in Fig. [2] The clean sys- 
tem with e = and hard wall boundary conditions has 
transmission T — 2 due to a non-chiral Majorana mode 
being present at each edge. Despite the fact that disor- 
der reduces this transmission, it only disappears after the 
bulk goes through a delocalization transition to become a 
trivial Anderson insulator at a sufficiently high disorder 
strength. 

To test the properties of the transmitting edges, we 
calculate the dependence of the transmission through the 
edges on the number of Kitaev chains N at a fixed e. In 
order to verify that it is the parity of the number of gap 
closings that determines whether an ensemble is an STI 
or not, we compare this result to the same calculation 
done for a mirror symmetric stack of chains in symmetry 
class BDI with alternating topological invariants Q — 
±2. Our results are shown in Fig. [3j We find that T ~ 
N^ 1 / 2 for the stack of Kitaev chains, as expected for a 
chain of randomly coupled Majorana fermions, or more 
generally for ID systems at the critical point [22l425j . 
Also in agreement with our expectations, we find that the 
edges of a BDI stack arc localized with T ~ exp(— cN). 

To test STIs in a dimension two higher than the di- 
mension of the original topological invariant, we consider 
a triangular lattice of Majorana fermions [2"rj] - f2"5] . which 
is a surface model of a 3D array of coupled Kitaev chains. 
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FIG. 3. Transmission as a function of number of Kitaev 
chains (blue, solid line) or of number of BDI chains with 
Q = ±2 (red, dashed line) in a stack. The transmission 
through the Kitaev chain stack scales as 1/yN, while trans- 
mission through the BDI stack is exponentially decaying. The 
inset shows the same curves in log-linear scale. In both cases, 
each point was obtained by averaging approximately 2 x 10 4 
disorder realizations at a fixed disorder strength e — 3t and 
chain length of 40 lattice sites. 



The tight binding Hamiltonian of this lattice is given by 

H = ^ity7t7j) (3) 
(y) 

with tij — —tji real- valued couplings between Majorana 
fermions. In a clean translationally invariant system with 
one Majorana fermion per unit cell and an arbitrary 
choice of the hopping signs, if the hoppings have equal 
magnitude, then the system has a reflection symmetry 
with respect to a plane passing through one of the hop- 
pings and a reflection antisymmetry (H changes sign) 
with respect to a perpendicular plane passing through 
any site. Since there is a single Fermi surface in the clean 
system, this means that any disorder that preserves the 
two mirror symmetries should make the translationally 
invariant triangular lattice of Majorana fermions a sur- 
face model of an STL 

We have calculated the dependence of the total trans- 
mission through a square-shaped region of a disordered 
Majorana lattice as a function of system size for different 
types of disorder. Our results are shown in Fig. |4j We 
have found that conductance scales positively with sys- 
tem size for any disorder strength both for log-normal, 
and uniformly distributed disorder, confirming our ex- 
pectations. Our results explain the thermal metal ob- 
served in Refs. [37J [25] • These two works numerically 
studied a transition between two triangular Majorana 
lattices with Chern numbers ±1 by changing the prob- 
ability distribution of hopping signs. They reported a 
metallic phase centered around the point where every 
hopping had a random sign. Since in that case the en- 
semble has a sufficiently high symmetry, and in particular 
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FIG. 4. Top left: Arrows show the directions in a triangu- 
lar Majorana lattice in which we choose the average hopping 
values to be positive. Bottom left: The setup in which we 
calculate transmission through disordered Majorana lattice. 
The leads are shown as grey rectangles. Right: Transmission 
through the Majorana lattice as a function of system size. 
The red curve corresponds to the case of hoppings tij uni- 
formly distributed in a range [— t, 5t\. The other three curves 
correspond to log- normal distributed tij, with zero log-mean 
and different values of the log- variance a. 



since it has the mirror symmetries described in the pre- 
vious paragraph, the metallic phase is a consequence of 
it being a surface model of an STL 

Conclusions and discussion. To conclude, we have uni- 
fied several known examples of systems protected from lo- 
calization by ensemble symmetries into the new class of 
statistical topological insulators. We were able to present 
a universal proof of why STIs avoid localization, by show- 
ing that the ensemble symmetry prevents them from hav- 
ing a definite value of a surface topological invariant. We 
have introduced a universal construction of STIs from 
any topological invariant using reflection symmetry, and 
were able to explain the thermal metal phase of Refs. 
[27l [28] as being a surface model of an STL 

Despite only having discussed non-interacting models, 
STIs should be protected from interactions, as long as 
the interactions do not introduce spontaneous symmetry 
breaking with long-range correlations. This holds since 
the identification of the topological invariant with a gap 
closing relies on the Green's function structure [20] . and 
is thus insensitive to interactions. For weak topologi- 
cal insulators, the robustness against interactions was 
reported in Ref. A natural extension of our ap- 

proach would include providing a more complete relation 
between ensemble symmetry groups and STIs. While 
we have focused on Z2 symmetries for simplicity, trans- 
lational symmetries must also be sufficient to construct 
STI for Z2 topological invariant, as is the case for the 
weak topological insulators, or also the triangular Majo- 
rana lattice. The STIs using fractionalized phases may 
provide a new way to make fractional topological insula- 
tors. The way the presence of several topological insula- 
tor and STI phases in the same symmetry class enriches 
the phase diagram of Anderson insulators provides an- 
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other interesting direction to study. Finally, we are as 
yet unable to solve the problem of efficiently evaluating 
the STI topological invariant on a general basis. It is suf- 
ficiently simple for several classes, e.g. weak TI, but be- 
comes hard for more complicated ensemble symmetries. 
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